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Abstract: We determine a next-to-leading order result for the shear channel thermal 
spectral function in SU(iV) Yang-Mills theory, working in the limit of vanishing external 
three- momentum. The result is subsequently applied to the evaluation of the corresponding 
imaginary time correlator, and its use in the context of sum rules is discussed. Our hope 
is that the calculation will eventually find use in the nonperturbative determination of the 
shear viscosity of the theory. 
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1 Introduction 

In recent years, a lot of effort has been put in the determination of various transport 
coefficients of the quark gluon plasma (QGP), characterizing its relaxation towards thermal 
equilibrium. This work has been largely motivated by the experimental heavy ion programs 
taking place at RHIC and the LHC [1, 2], which have highlighted the importance of not 
only understanding the bulk equilibrium properties of the produced matter but also its 
hydrodynamic behavior. One particularly prominent observation in this direction has been 
the discovery of the sizable effect that the small, yet finite value of the shear viscosity r\ 
has on the flow properties of the plasma [3, 4]. 

Unfortunately, the first principles determination of transport coefficients in an inter- 
acting quantum field theory is a notoriously hard problem. In the case of heavy ion physics, 
a further complication arises from the fact that the matter produced at RHIC has been 
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observed to behave more like a strongly coupled liquid than a gas of weakly interacting 
quasiparticles, clearly invalidating its description in terms of perturbative methods (which 
furthermore are extremely difficult in this context, cf. e.g. [5, 6]). At the same time, a 
reliable nonperturbative determination of transport coefficients, in principle available from 
the zero frequency limit of the corresponding (Minkowskian) spectral functions, has proven 
extremely challenging, as the lattice approach is restricted to the Euclidean formulation of 
the theory [7-9]. As a result, a considerable amount of attention has shifted towards the 
gauge/gravity duality, which has in particular led to the famous claim of the QGP being 
a nearly 'ideal' fluid, with a shear viscosity to entropy ratio close to the conjectured limit 
v/ s — V^ 71 ") [10] (see also [11, 12] for some interesting recent developments). 

Clearly, it would be desirable to be able to determine parameters such as the shear 
viscosity of the QGP using nonperturbative tools but staying within the physical theory. To 
this end, an ambitious program was initiated in [13], where the authors tested a numerical 
recipe for the analytic continuation of Euclidean lattice data to Minkowskian signature. 
This method in practice amounts to the inversion of an integral relation between the 
spectral function p(oj) and the corresponding imaginary time correlator G(t), 



constituting an in principle well defined but numerically extremely challenging task. While 
early results have been encouraging (see also [14, 15]), they have also highlighted the 
necessity of subtracting short-distance divergences from the Euclidean correlators and in 
general obtaining as much analytic information on their behavior as possible. This in- 
cludes not only determining the temperature independent ultraviolet (UV) singularities, 
available from T = perturbation theory (see e.g. [16]), but also the leading temperature 
dependent terms. The latter quantities are most conveniently obtained from the behavior 
of the corresponding spectral functions at uj ~ irT, computable using finite-temperature 
perturbation theory. 

For the bulk channel of SU(iV) Yang-Mills theory, including operators such as the 
trace of the energy momentum tensor, the task of perturbatively determining the spectral 
function was undertaken in [17]. There, the quantity was evaluated to next-to-leading 
order (NLO) in the limit of vanishing external three-momentum, building on the earlier 
work of [18, 19]. In the paper at hand, our goal is to extend this result to the perhaps even 
more interesting but also technically more complicated shear channel, which necessitates 
the generalization of the methods of [17] to a number of new master integrals, involving 
components of the loop three-momenta in the numerator. As a preparation for this, the 
UV limit of the shear correlator was analyzed in terms of an Operator Product Expansion 
(OPE) in [20], leading to results compatible with known sum rules [21, 22] as well as the 
arguments of [23]. Our hope is that the full NLO thermal spectral function will eventually 
find use in the extraction of the shear viscosity, 
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from lattice data. In addition, our results can immediately be used to evaluate the Eu- 
clidean imaginary time correlator, a quantity that has been addressed in several recent 
lattice and AdS/CFT calculations in both the shear and bulk channels (see e.g. [24-29], 
some of which also address spatial correlators). 

The outline of our paper is as follows. In Section 2, we present our setup and introduce 
the necessary notations, while Section 3 contains an overview of how our computation was 
performed. Section 4 then lays out our result for the spectral function and discusses its 
implications for the Euclidean imaginary time correlator and sum rules, after which we 
draw our conclusions in Section 5. Most details of our calculation, including the definitions 
of all master integrals, can be found from appendices A and B. 



2 Setup 

The setting of our calculation is very similar to that introduced in detail in [17, 20], and 
for simplicity we adopt our notation from these references. We work within pure SU(iV) 
Yang-Mills theory at a nonzero temperature T, defined by the dimensionally regularized 
Euclidean action 



1_„ _„ f 1 



with f3 = 1/T and 



At I d^x^F^ ee / \F; U F%, (2.1) 



F; v ee d^Al - d v Al + g B f abc AlAl . (2.2) 



It is straightforward to see from here that the energy-momentum tensor of the theory takes 
the form 

1 

I' 
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the (connected) correlation functions of which we set out to compute 

G^ a p{x) ee (T^(x) T a p(0)) c . (2.4) 
In momentum space, the corresponding Euclidean correlator takes the form 

G^ aP (P) ee [ e- ip - x G^ aP {x) , (2.5) 



where the Fourier transform is taken in D = 4 — 2e dimensions and from which the desired 
spectral functions can be obtained, 

P^M") = Im G»vM P ) „ , ^. n+ln , • (2-6) 

As discussed in some length in [20], to evaluate the momentum space shear correlator 
^12,12 (P) within dimensional regularization, it is convenient to introduce the projector 
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Figure 1. The one- and two- loop Feynman diagrams contributing to the NLO shear spectral 
function in SU(7V) Yang-Mills theory, with the curly line corresponding to the gluon field. The grey 
blob in (v) denotes the one-loop gluon self energy. 

where Pj^ u (P) is a symmetric transverse projector orthogonal to the four- vectors P and 
U = (1,0), satisfying Pq = = P^, P? = 5{j — (for convenience, we keep p nonzero 
here). Using this definition and choosing the spatial momentum p to point in the xd-x 
direction, we obtain 

X^af, G ma p(P) = —D(D - 2){D - 3) Gi 2)12 (P) . (2.8) 

This suggests that if we define the momentum space shear channel correlator through 

G V (P) = ZX^ap G^ U)a p(P) , (2.9) 

then in the e — > limit the corresponding (p = 0) spectral function p^oS) satisfies 

Pn(u) = -16pi2,i2 (w). (2.10) 

Finally, we note that our calculation will be carried out within the MS renormalization 
scheme. The sum-integral measure we use reads 



where A and A stand for the renormalization scales in the MS and MS schemes, respectively. 
3 Calculations 

We begin the determination of the NLO spectral function p n [u) from the Euclidean mo- 
mentum space correlator G ri (P), which can be seen to obtain contributions from the one- 
and two-loop graphs displayed in fig. 1. Taking advantage of eq. (3.2) of [20] and drop- 
ping contributions that depend on the external momentum in a trivial way and thus don't 
contribute to the spectral function, we see that this quantity can be written in the form 

= PjS^ + ^c{px b H + Px d M + px f {u) + pxjyi) + p X; (w) + Px . (w)} , (3.1) 
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where cIa = N? — 1 and we have denoted 



Pj h (" 



Px d i u 

Px^ 
Pt x {u 

Px.(u 
Px. (w 



D{D - 2)(D - 3) 
8 



Pj (uj) - D(D - 3)pj 



(D-2)(D-3)p j2 (co) , 
D(D-2){D-3) 



2p xl (u) + p x2 (uj) 



= D(D -2){D- 3) \o z i (u) + 3p x2 {u) ) + 
Dp x0 (lo) + 8p X 3 (w) 

d d 

3D{D -2)(D -3) 
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P x i{uj) , 

* ( \ a f \ f \\ D{D-2) 
Ap x o (u) + 4 % i (w) - ^3 (w) j + -^— 2 '- 



2p xl (u) + 2p X 2 (uj) + p x? (u/ 



p x o (w) + -D(£> - 3) 2/o x i (w) + p X 2 (w) 



(3.2) 
(3.3) 

(3.4) 
(3.5) 

(3.6) 
(3.7) 



+ 



12 - 16D + 3D 2 



2 Px? (uj) + p xA {oo) - D{D - 6) p x5 (uj) + 2p xp (u) . (3. 



The definitions of the various 'master spectral functions' appearing here can be read off 
from appendix A.l; note that our notation for them differs somewhat from that of [20]. 
The primary computational task of the present paper will be to derive explicit results 
for these quantities, evaluating the divergent terms analytically and the most complicated 
finite parts numerically. 

The general techniques required in the computation of the above masters were devel- 
oped and rather thoroughly explained in [17]; in fact, many of the results needed in the 
present calculation can be directly read off from this reference. The main complication 
involved in the evaluation of the masters only appearing in the shear channel is related to 
how one should deal with the projection operators Pt{Q) = Q^QuP^ u (P), which appear 
frequently in the definitions of appendix A.l. Here, a crucial simplification appears due to 
setting the external three- momentum p to zero, which allows us to average the results over 
its direction, applying the simple relations 



1 



D - 1 
1 



Sij , 



(3.9) 
(3.10) 



(PiPjPfcPz)p = D 2 _i (Sijhi + $ik$ji + SuSjk) . 
This leads to the replacement rules (with Q and R standing for Z?-dimensional loop mo 
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menta here) 



p T {Q) §^ 2 , pKq) ^r^ 4 ' 

p r (Q)^(i?) D2 ~ 2 2 ^' 2 ^ + • r) 2 , (3.H) 



which will prove extremely useful in the following. 
Next, we introduce the notation 



Pjj?( u ) = / /J?. P2$H= / > ( 3 - 12 ) 

Jq J q,r 

in which the one- and two-loop functions fjo and f%o, corresponding to the integrals en- 
countered in the bulk case, can be read off from [17] and are for completeness reproduced 
in appendix A. 2. Using eq. (3.11) then allows us to immediately write down similar ex- 
pressions for several of our new masters, 



, D-2q 2 _ D(D - 2) g 4 

hi --flTT^ b ° ' JJg- D 2_ 1 ^Jjg> 

_ D-2 q 2 D-2 r 2 

D-2 q 2 

L d D- \ u 1J - L >- ' J1 

lh.l) 2),' , _, 



/z3 - ~ — 7772/zS ' /z2 - -77 — 7772/15 ' 



D - 


Ioj 2 


D - 


2 r 2 


D - 


Ioj 2 


D - 


2q 2 


D - 


lco 2 


D - 


2r 2 



D - 2 q 2 D-2E 2 r -\ 2 
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D(D - 2) g 4 - 2) - A 2 ) 2 
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1 ^ D 2 -2D- 2 2 2 2 2 \ 
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3.13) 
3.14) 
3.15) 
3.16) 
3.17) 
3.18) 
3.19) 
3.20) 
3.21) 



where A is a fictitious mass parameter that was introduced in [17] to regulate IR divergences 
present at intermediate stages of the computation (but which will all cancel in the end, 
allowing for the A — > limit to be taken). 

The remaining master integrals — X^ , X 4 , Xf , I®, Xj, X*, X 2 , if — differ from the ones 
evaluated in [17] not only by the presence of projection operators but also because they 
each contain a squared propagator. To deal with the complications induced by this, we 
have found it most convenient to introduce a mass parameter m in this special propagator 
(chosen to carry the loop momentum R), evaluate most parts of the integrals with m^O, 
and in the end make use of the identity 

w = - lim < -A?^r^ — o \ ■ (3-22) 
i? 4 m->0 I dm 2 R 2 + m 2 I V ' 
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Including this mass into the definition of 2^ , we then quickly arrive at the following results 
for the h type masters: 



IT* 




(3.24) 

H = ^ T 1™ 4 ihiTTih* ) > (3-25) 



^ = & { ^ r ^ - 1 v (q - r) ^ + ^r (q • (q - r))2 ) S I • (3 - 27) 

An important simplification in the evaluation of these integrals comes from the fact that 
after the introduction of the parameter m, the IR regulator A is no longer needed. A 
nonzero value of m is seen to regulate all IR divergences in the above expressions and lead 
to a finite result when the h type masters are added together and the m — > limit taken. 

Moving finally on to the i type integrals, the introduction of the parameter m in the 
i?-propagator of IP leads to the decomposition 



1° 



^ 4(Q-P) 2 + m 2 P 2 ^ m 2 

\ R Q 2 [R 2 + m 2 }(Q - R) 2 (R - P) 2 + % 



q,r Q 2 [R 2 + m 2 ](Q - R) 2 (R - P) 2 %^Q 2 [R 2 + m 2 ]{R-P) 



2 



p2 p p2 



^2 



^ R Q 2 (Q-R) 2 (R-P) 2 r Q , R Q 2 [R 2 + m 2 ](Q-R) 
+ t Q , R Q 2 [R 2 + m 2 } ~ t Q , R Q 2 ^ 2 ' (3 - 28) 

where only the first three terms give nonzero contributions to the spectral function. A 
straightforward calculation now leads to the result given in eq. (A. 39) for fji, while the 
other i type integrals become 

D-2 ,. [ d 2 ,1 

With the above relations, the evaluation of the master integrals — and thus the entire 
NLO spectral function — becomes an in principle straightforward, though technically quite 
demanding exercise in integration. In order not to burden the reader extensively, we leave 
the details of the practical implementation of these integrals to appendix B, and instead 
move directly to the results. 
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Figure 2. The behavior of the function </>^(w/T) on a linear and a logarithmic scale, multiplied 
by (uj/T) 3 in the former case. 



4 Results 

4.1 Spectral function: Limits and numerical evaluation 

Results for each of the functions p x (u>), appearing in eq. (3.1), are given in appendix B. 
With the exception of a few very simple cases, they are composed of two parts: An analytic 
piece containing the possible UV divergences and logarithms of the renormalization scale, 
and a numerically evaluated finite part denoted by px and consisting of (often rather 
complicated) one- and two-dimensional integrals. Collecting these results together and 
inserting them to eq. (3.1), we see that the entire NLO shear spectral function becomes 

^^ + ^-^(^H}, ,4,) 

where n x = l/(e@ x — 1) and we have defined the dimensionless function (cf. eqs. (4.1)-(4.2) 
of [17]) 

4(u,/T) = " h^/T) + Px^/T) - px<(u/T) - px^/T) . (4.2) 

This function is plotted in fig. 2, and the corresponding lengthy integral expression repro- 
duced in the Mathematica file shearresults.nb, available at [30]. 

Before proceeding to plot the above result, let us first perform a number of consistency 
checks on it. The first thing to note is the automatic cancelation of the 1/e divergences in 
eq. (4.1), as well as the agreement of the constant terms with the T = results of [16, 20] 
(note in particular that no renormalization is required at this order). Inspecting next the 
subleading terms at large co/T, we note that the functions px x (u/T) can be analytically 
expanded in inverse powers of T 2 /lo 2 . This results in the expressions 

Px>/T) = ^f + + 0(T^ 6 ), (4.3) 



3u 2 9w 4 
38tt 2 T 2 392vr 4 T 4 
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^/ T > = ^ + -mjr +00-*). <«> 

which — when inserted into eq. (4.2) — exactly cancel each other to the order indicated. 
In fact, it can be shown that in this limit of oj — > oo the function (p^(oj) behaves as 

#.( W /T) = i^- + 0(T 8 /^ 8 ) , (4.7) 

which confirms the expectations of [23], and in addition constitutes a highly nontrivial 
cross check of our result. 

In the opposite limit of small oj/T, our result is dominated by two integrals, px h (oj/T) 
and px. (oj/T), which have the IR expansions 

p Xh (u/T) = + 0(T> 2 ) , (4.8) 

3^7r 2 T 3 

PtHT) = -^3- + ^(r 2 /^ 2 ) • (4.9) 
This in particular implies that the full function <j>j<(ui) behaves at small oj as 

= 7 3 +O(r 2 /o; 2 ), (4.10) 

and that the NLO shear spectral function approaches a (positive) constant of 0(g 2 ) in the 
w — > limit. This is in clear contrast with the corresponding limit of the bulk spectral 
function considered in [17], and obviously does not represent the physical IR behavior of 
the quantity. 1 For oj = O(gT), a naive perturbative expansion of the type of eq. (4.1) 
namely breaks down, and a resummation of the soft contributions to the spectral function, 
available through the Hard Thermal Loop (HTL) effective theory, must be performed. For 
the bulk case, this lengthy exercise was carried out in [17], but in our current work, we have 
decided to refrain from it. The reason for this is that our result is in any case expected 
to be trustworthy (and find its applications) in the regime oj > T, where the contribution 
of the extra HTL terms is negligible (cf. fig. 4 of [17]). In addition, as will be discussed 
in the following subsection, the contribution of the perturbative IR sector to quantities 
such as the imaginary time shear correlator is expected to be numerically rather strongly 
suppressed. 

With the above reservation, we now proceed to numerically evaluate the spectral func- 
tion by inserting the function (^(oj/T) into eq. (4.1). This leads to the behavior displayed 
in fig. 3, where we have used a one-loop running coupling, set the number of colors iV c to 
3, and varied the renormalization scale A by a factor of two around the 'EQCD' value [31] 

In h& = _ 7E _ 1 . (4.11) 
4ttT ' 22 v ' 



1 It is nevertheless interesting to compare this result to the findings of ref. [11], where the authors report 
that the physical w — > limit of the spectral function is a constant of 0(g 8 ). 
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Figure 3. The behavior of the absolute value of the shear spectral function, |p^(w)|/(4d^T 4 ), for 
T = 3T C (corresponding to 3.75A^g). The two blue curves stand for the NLO result evaluated with 
A = 0.5A opt and A = 2A opt , while the dashed red curve shows the leading order (LO) result. The 
spikes in the NLO results correspond to a change of sign in the quantity from positive to negative 
(with increasing oj). 

As there are no explicit logarithms of the renormalization scale in eq. (4.1), an optimization 
scheme such as the Fastest Apparent Convergence is clearly not available for us. At large 
values of lj, a choice of the form A ~ oj might seem more natural than the above, as it would 
lead to a logarithmically vanishing coupling; as seen from fig. 3, even a scale proportional 
to T, however, leads to a fast convergence of the result towards the free theory limit. 

4.2 Sum rule and imaginary time correlator 

As in the case of the bulk channel [17], there are two natural immediate applications of 
the shear spectral function: The verification of the sum rule of eq. (4.19) of [20] (see also 
[21] and [22]), and the determination of the imaginary time shear correlator, cf. eq. (1.1). 
This time there are, however, two problems preventing a straightforward evaluation of the 
integrals appearing in these relations. First, in the oj — > limit, the shear spectral function 
approaches a constant, rendering both of these integrals logarithmically IR divergent. And 
second, somewhat connected to the first issue, already at LO the spectral function is known 
to obtain contributions of the form oj5(oj), which do not show up in our calculation but 
nevertheless give nonzero contributions to integrals of p v (oj)/oj. At two loop order, some 
or these terms are in addition proportional to IR divergent momentum integrals, further 
complicating their treatment. 

Our approach to resolving the above issues is somewhat different in the two cases. 
For the sum rule, we expect that a successful verification of the identity will be possible 
only after extensive further work on the oj — > limit of the spectral function, involving an 
HTL type resummation. As we do not wish to embark on such a calculation here, we can 
verify the sum rule only for those master integrals, for which the IR issues do not show 
up (note that in the form of eq. (4.19) of [20], the sum rule should indeed hold for each of 
the master spectral functions separately). Fortunately, this class of masters includes most 
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Figure 4. The imaginary time correlator G5j cf (r), denned by eq. (4.12) and displayed on a loga- 
rithmic scale. Just as in fig. 3, altogether three curves are displayed here, but as is apparent form 
the figure, the two blue NLO curves practically overlap with the dashed red (LO) one. 

of the integrals considered in this work, in particular all of the j types and most h's as 
well. In all of these cases, we have successfully verified the sum rule to the accuracy of our 
numerical integration, which serves as a powerful check of our results. 

Moving next on to the imaginary time correlator, we again expect to be able to do 
the full integral in eq. (1.1) only after a proper handling of the IR limit of the spectral 
function. However, our expectation based on fig. 3 as well as fig. 4 of [17] is that the 
physical contribution of the very softest frequencies to this quantity will most likely be 
numerically fairly small. Neglecting the constant shift of the correlator produced by the 
unknown uj5(u) contributions, we may thus define a deformed correlation function 



where loq ~ 0.6T stands for the frequency, at which the NLO spectral function changes its 
sign. This is in fact a very natural choice, as below this value the LO spectral function is 
seen to obtain a relative correction of more than 100%. 

The 'deformed' imaginary time correlator is plotted in fig. 4, which displays an excep- 
tionally good agreement between the LO and NLO results. This behavior can be attributed 
on one hand to the small relative size of the NLO correction in the T = spectral function, 
cf. eq. (4.1), and on the other hand to the vanishing of the leading temperature dependent 
corrections to the quantity at large to. It would be very interesting to quantitatively com- 
pare this behavior to existing lattice and AdS predictions (see e.g. [7, 29]), which in fact is 
a question that will be addressed in the near future [32] . 

5 Conclusions 

In the paper at hand, we have performed a NLO calculation of the thermal spectral function 
in the shear channel of SU(iV) Yang-Mills theory. The main results of the paper were 




(4.12) 
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presented in sec. 4, with the most important finding being the behavior of the function 
4>j,(u>/T) that encodes the temperature dependence of the spectral function. In accordance 
with the general arguments of [23] as well as the OPE results of [20], we found that 
this function vanishes in the UV limit like 1/w 6 , with a proportionality constant given 
in eq. (4.10). This result was subsequently used to determine the imaginary time shear 
correlator, in which the NLO term was seen to be highly suppressed. 

On the computational level, the derivation of the NLO spectral function necessitated 
an extensive generalization of the methods developed for the bulk channel calculation of 
[17]. Apart from the the results displayed in the previous section, a very important part of 
our paper can indeed be found from sec. 3 as well as appendix B, which detail the techniques 
we used in our work. Most importantly, these sections contain a detailed explanation of 
how one can efficiently deal with components of three-momenta as well as squared bosonic 
propagators in the integrands of two- loop master spectral functions. 

In the future, we hope that our results will find important use in the eventual nonper- 
turbative determination of the shear viscosity of Yang-Mills theory, following the rationale 
of [14]. To this end, our finding concerning the small magnitude of the finite temperature 
corrections to the spectral function — and in particular to the imaginary time correlator 
- is clearly intriguing. It in particular highlights the relative importance of the sub- 
leading corrections to the T = spectral function, which incidentally have recently been 
determined to an impressively high perturbative order [16]. 

Finally, we note that amongst the directions, to which our present calculation should 
be continued, two stand out quite clearly. First, for the sake of completeness, it would be 
important to be able to properly determine the IR limit of the NLO shear spectral function, 
so that all contributions of this order to the sum rule and imaginary time correlator could 
be accounted for. This would imply both an HTL type resummation at frequencies uj ~ gT 
and a successful treatment of the uj8{uS) contributions to the spectral function. At the same 
time, it would also be interesting to generalize the present calculation to the presence of 
(massless) fermions, extending the theory to full QCD. In this case, the result is expected 
to contain T-dependent terms of orders l/ta 2 and 1/w 4 [23], so a natural first step might 
be to generalize the OPE results of [18, 20] to the presence of quarks. 2 
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A Master integrals 

In this first appendix, we first list all master integrals appearing in our calculation, and 
then recall from [17] the forms of the integrand functions fj , f x0 , defined in eq. (3.12). 

A.l Definitions 

The master integrals appearing in our result (3.1) for the shear spectral function read 
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In these expressions, we have defined P T {Q) = Q^QuP^ u {P) = q 2 - (q ■ p) 2 , where P^ U {P) 
is the transverse projection operator discussed in section 2. 

A. 2 Integrands 

In section 3, we reported a number of relations between the integrands f%n of different 
master spectral functions, relating them to cases encountered in [17] (albeit with a different 
IR regulator). Using a notation, where E r = \/ q 2 + m 2 , E qr = |q — r|, and n qr = riE qr , 
the latter read: 
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B Evaluation of the masters 



In this appendix, we provide details on how the master spectral functions of eq. (3.1) were 
evaluated. We do this case by case, starting from pj and ending with p x . Note that in 
this process, we will as a rule discard all contributions that would lead to contact terms 
proportional to 5{ui) or derivatives thereof. 
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Summing up the different parts contributing to pj , we finally reach the compact result 
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this immediately gives us 
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For p xX , we on the other hand easily obtain 
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Finally, the corresponding expressions for /) t2 read 
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Starting with p.^, we follow section B.6 of [17] and obtain 
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For p Tl , we on the other hand first keep m nonzero, obtaining 
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With p l3 , we finally get 
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which leads to the desired result 
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Next, we move on to master integrals of type h, which represent our first truly two-loop 
topology. We begin from a detailed treatment of p x o , which has been considered in section 

h 

B.10 of [17], but is now generalized to the presence of the mass parameter m. As discussed 
above, this parameter will not only enable us to derive results for other master integrals 
containing squared propagators, but its nonzero value in addition serves as an infrared (IR) 
regulator, implying that we may set the A parameter of [17] to zero. For brevity, we will 
write down an explicit final result only for p x o, while for the other integrals, we simply 

h 

indicate, how the calculation proceeds. 

The evaluation of p r0 begins from eq. (A. 37), obtained by performing the corresponding 
Matsubara sum with nonzero m. From here on, we follow the treatment of [17] and divide 
the terms appearing in this expression to three categories, dubbed 'factorized powerlike' 
(fz,p), 'factorized exponential' (fz,e), and 'phase space' (ps) integrals. For a discussion of 
the physical interpretation of these contributions, see appendix A.l of [17]. 
The (fz,p) contribution corresponds to the UV divergent expression 
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where again E qr = |q — r| and the r integral can be easily performed using 
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The remaining integral over q on the other hand reduces to 
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In the last stage, we have set ng r + n r = 2n^ , owing to the identity 
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For the UV-finite (fz,e) part, corresponding to the remaining terms on the first six 
lines of eq. (A. 37), we perform a change of integration variables according to 
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where P stands for a principal value integral and we have used the fact that the integrals 
on the two first rows of eq. (B.28) are clearly identical upon the redefinition q -H- E qr . 

The (ps) part of the integral, corresponding to the last four rows of eq. (A. 37), is 
technically the most complicated one. For it, we begin by writing the integration measure 
in the form 
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where E qr are defined as above. Following then the steps laid out in section A. 4 of [17], we 
reach (after some labor) the result 
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where we have made use of the symmetries of the integrand to write 
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in the first two parts of the expression. In addition to this, we have for further reference 
introduced the function F(x, y, z) in all of the integrands, which in the present case obtains 
the value F x o(x, y, z) = 1 . It is straightforward to see that eq. (B.30) reduces to eq. (B.42) 
of ref. [17], when the limit m — > is taken. 

Next, we must perform the integrals in the above three parts (i)-(iii). In each case, we 

4 

separate the result (with the overall normalization factor p^w (1 + 2n«) left out) to three 
subcontributions based on the distribution functions they contain, denoting by A those 
parts either containing one n q + r but no other re's or unity; by B those with distribution 
functions independent of r; and by C the remaining piece, proportional to a negative power 
of n r . In some cases, we will also find it convenient to perform a change of variables from 
q and r to q + r = x and q — r = y. 

(i) Taking advantage of the symmetry of the integral, it is easy to see that the A contri- 
bution can be written in the form 
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in which we have at the last stage renamed x = q and where the 'vacuum' contribution 
is clearly analytically integrable. For the B and C parts, we use the symmetry of the 
integrand to restrict the integration region to the q > r part, ending up with 
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The latter of these expressions must be evaluated as a two-dimensional numerical integral, 
(ii) For these integrals, the exact same steps as above lead to the results 
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Collecting all of the above results together and organizing the terms somewhat, we 
finally obtain a lengthy expression for the function p x0 (w) 
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Unlike in the (fz,p), (fz,e) and (ps) parts separately, taking the m — > limit of this 
expression leads to a finite result, which in fact is seen to agree with that derived in 
ref. [17]. We will nevertheless keep m nonzero for the time being, as this will turn out 
useful in the following. 
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As discussed in sec. 3, the integrand of p Tl is related to that of p T0 through 
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from which it is easy to derive the (fz,p) result 
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For the (fz,e) part, we similarly obtain 



(fz,e) / N 



UJ 



2(4tt)3 C 1 + "*•+"■■ 



+ 





oo 



d<? / dE qr n q q , 



1 1 

+ 



1 



A 



10 



A 



01 



qr 



r+E„ 



dq n qr q 



\r-E qr \ 



A 



+ 



00 



A, 



+ 



01 



A 



+ 



An 
1 



10 



A 



n 



(B.42) 



.(B.43) 



Proceeding finally to the (ps) part, we recall that in the derivation of the equivalent of our 
eq. (B.30) in [17], a change of variables of the type 
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was carried out. It is then straightforward to see that the (ps) contribution to X\ takes the 
form of eq. (B.30) with 
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All integrals encountered in evaluating this expression can be dealt with using the methods 
discussed above. 
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D-2r 



i 2 



2 H ' 



D-Iuj 

from which one straightforwardly obtains the relations 
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as well as 
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Again, all integrals encountered are carried out with methods used in the previous cases. 

With 1^ , we encounter the first master integral with a squared propagator. To this end, 
we take a mass derivative of p x2 (oj), obtaining 
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As we have above computed p x i{oj) while keeping m nonzero, taking the derivative is in 
principle a very straightforward task. The only somewhat problematic issue is related to 
the IR divergences that appear in individual parts of the integral once one proceeds to set 
m — > in the end. Our strategy with them is to identify those parts of the integrals that 
diverge in this limit, and subseqently add and subtract from the integrands terms that are 
analytically computable (with finite m), yet render the original integral convergent. The 
divergent terms are expected to cancel against each other once we assemble the result for 
the full master spectral function, while the remaining finite parts are dealt with using the 
methods described above. 

In the one-dimensional integrals, originating from the A and B parts defined above, 
the divergences in the m — > limit may in principle appear due to one of two reasons: A 
pole at q = ^ , coming from a factor l/(q — | + ^ ), or an explicit lnm 2 term originating 
from the analytic r integral. For , the latter does not occur, while the former can be 
dealt with by writing 
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Here, the integral containing the curly brackets is seen to be finite, while the divergent 
term on the last row can be evaluated analytically, producing 
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Turning our attention next to the two-dimensional part, we observe that the only 
divergence originates from the r = limit of This can be regulated by subtracting 

2 

from the integrand a term of the type a/{r + tt-), where a is a suitably chosen residue. 



This amounts to writing in the somewhat complicated form 
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where the two terms on the latter row combine to a finite integral. 

Apart from the above terms, all other parts of p l3 are IR finite and computable with 
standard methods. 



This time, we begin with the integrand relation 
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using which we may define a non-differentiated two argument version of the master integral 
p l4 (uj,m), 
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as well as 



TP ( \ 



(B.57) 



where the function F x a{x, y, z) refers to the non-differentiated version X 4 as well. The 

h 

differentiation of these expressions with respect to m and the subsequent subtraction of 
divergences proceeds exactly as explained in the previous cases. 
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For 1^, we again use the definition 



D(D-V llm d H (B58) 
to define the non-differentiated version of the integral, p T5 (a;,m), through 

For this function, we easily obtain the results 
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The master integral I® has a rather complicated structure in momentum space. It is related 
to X° via the relation 
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which we again write in the form 
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From here, we straightforwardly obtain the (fz,p) result 
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but to evaluate the other parts requires a little more care. 

To deal with the (fz,e) and (ps) contributions, we first note that the full master can 
be written in the form 
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For the (fz,e) parts of these integrals, we then obtain 
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while the (ps) contributions are available via 
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The evaluation of the individual integrals is again a relatively straightforward exercise, after 
which the final result for the master is available by simply adding the different contributions 
together. 
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Figure 5. The behavior of the absolute value of pi (oj/T), displayed on a logarithmic scale. The 
spike in the curve corresponds to the sign of the function changing from negative to positive with 
increasing uj. 
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This time, our starting point is 
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which leads to the rather convenient result 
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For this function, the (fz,p) contribution is easily evaluated and reads 
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while the (fz,e) and (ps) contributions are immediately available using previous results. 
Final result 

Collecting all of the above pieces, proceeding to the m — > limit, and performing the 
remaining convergent integrals numerically, we arrive at the result 
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Here, px h (uj/T) is a dimensionless quantity, whose behavior as a function of frequency we 
display in fig. 5. 

B.5 p X{ 

Moving on to X f , the only master integral needed is p Tl (w), related to the already evaluated 

p (w) (see section B.8 of ref. [17]) through 
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Noting that this integral only contains a (ps) part, we quickly obtain the convergent result 
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where we have made use of the change of variables introduced in eq. (B.44). Utilizing 
finally the e = relation 
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this provides us with the final result for p x (u), 
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Here, the constant term corresponds to the T = limit of the result, and the function 
Px f (w/T) is displayed in fig. 6. 

B.6 p x 

With the i type integrals, we start from the i', which does not appear in the expression for 
p x , but is required in the evaluation of p x? . 
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Proceeding as before, we begin by reading off the (fz,p) contribution from eq. (A. 38), 
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Similarly, the (fz,e) part gives 
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while the (ps) contribution can be related to the corresponding part of the p x ± integral 



leading us to the final result for the spectral function. 
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With if, it is a straightforward exercise to show that the integrand of the spectral function 
is related to that of p x (uj) through 
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which in combination with our previous results complete the evaluation of the integral. 
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where we have made use of the shifts of variables described in eq. (B.44). 
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With the results for X} at hand, it is very straightforward to to obtain the corresponding 
expressions for I 2 , utilizing the relation 
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A simple differentiation now results in the expressions 
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Final result 



Having evaluated all of the different pieces contributing to p x , we are ready to write down 
our final result for the quantity. Proceeding to the m — > limit, we see that all IR divergent 
terms cancel, leaving the result 

lt . . u A , . f 1 1 A 2 23 2tt 2 T 2 , , . 1 . 
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Figure 7. The behavior of the function pi. (uj/T). The spike corresponds to the sign of the function 
changing from negative to positive. 



where the dimensionless function p%. is displayed in fig. 7. 
B.7 p Xj 

The j type integrals represent the most complicated topology encountered in our calcula- 
tion. A minor simplification occurs, though, due to the fact that no squared propagators 
appear in these integrals, implying that instead of introducing the mass parameter m as 
above, we may use the IR regulator A introduced in [17]. As this parameter appears as 
the mass of the Q — R propagator, it doesn't break the Q R symmetry of the integral, 
which turns out to simplify its evaluation somewhat. 

We begin the calculation by briefly recalling how the simplest case p r0 was dealt with 

j 

in [17]; the formulae introduced here will turn out useful also in the other cases considered. 

Pip 

j 

Recalling that the integral p x o is UV convergent, we may set D = 4 from the beginning. 

j 

The (fz,p), (fz,e) and (ps) parts of the integral can then be directly read off from [17], with 
the results taking the form 
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where we have now denoted E qr = \/(q — r) 2 + A 2 . Note also that our notation for the 
integration variables in the (fz,e) part is somewhat different than that of [17]. 

The evaluation of the above integrals is explained in great detail in appendix A of 
[17]. For the (fz,p) and (fz,e) parts, the expressions can (after some labor) be reduced to 
one-dimensional integrals over elementary functions, whereas the (ps) part leads to both 
one- and two-dimensional integrals. The parameter A can be set to zero at the very end, 
upon an explicit cancelation of the IR divergent terms. 



In the evaluation of p x i(u;), our starting point is the simple relation 
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Unfortunately, the altered UV behavior of the integrand has the effect that unlike with 

p x o(uj), we cannot set D = 4 yet, but rather have to first subtract an analytically evaluat- 
j 

able UV divergent part from the rest. To this end, we first consider the (fz,p) contribution, 
where this subtraction must be implemented, obtaining 
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The most convenient way of evaluating the new integral involves the use of the identity 
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where we have taken advantage of the fact that p = 0. Using this relation, we obtain 
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where the UV divergence is now contained in a simple one-loop T = integral. 

To evaluate the above integral, we add and subtract from its integrand the leading 
terms of a large-Q 2 expansion, writing 
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where m r is an in principle arbitrary regulatory mass. After this subtraction, the integral 
inside the curly brackets becomes finite and can be carried out with the same methods 

we used when dealing with p^' p \u), cf. [17]. The subtraction term, on the other hand, 

j 

is evaluated analytically. Here, a potential problem originates from a pinch singularity at 

2 2 

q = - ~^™ r , which can however be moved to the limit q = by the choice m r = u, where it 
is protected by the integration measure. After a straightforward exercise, these steps lead 
us to the result 
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With the (fz,e) and (ps) parts, no divergences appear, and we may proceed as with 



T?. This quickly leads us to the results 
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With p X 2 j we start from the relation 
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which immediately produces the (fz,p) contribution 
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Of the three terms on the last line here, only the middle one is both UV divergent and 
new. For it, we again write 
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where the UV subtraction term clearly vanishes in the limit of our interest, 
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The evaluation of the (fz,e) and (ps) parts of the integral proceeds as before. After 
some work, we obtain the results 
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which can again be straightforwardly evaluated with the methods introduced in [17]. 
Starting with 
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Here, the new UV divergent integral reads 
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where the UV part gives 
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and the integrals inside the curly brackets are again finite and can be treated with numerical 
methods [17]. 

Finally, the remaining contributions to this master obtain the forms 
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which we deal with as before. 
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Next, consider p xB . For this integral, the defining relation reads 
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Here, the first UV divergent integral was evaluated in connection with Zj 1 , while the second 
leads to the expression 
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Finally, the (fz,e) and (ps) contributions to the integral read 
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The evaluation of p^u), denned by 

1 / D 2 — 2D — 2 2 \ 

/zf = ^ ( ^2_! ^ " A ') + WZT^ ■ (q " r )) 2 J /*? , (B.139) 

is somewhat simplified by separating from it two known contributions according to 
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This quickly leads to the (fz,p) result 
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where the second term on the last row is furthermore UV convergent. The first term, on 
the other hand, produces 
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Once again, it is a simple exercise to write down the two last parts of our integral, 
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The reason we have left the evaluation of p x ±{oS) as the last one becomes obvious when one 
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writes down its definition 
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in terms of the variables q and r. For the (fz,p) part, we namely obtain 
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where each of the UV divergent integrals in the one-loop expression inside the curly backets 
has been evaluated above, cf. eqs. (B.113), (B.120), (B.127), (B.134), and (B.142). 

With the (fz,e) and (ps) contributions, the calculation proceeds as in the previous 
cases, producing 
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Figure 8. The behavior of the absolute value of the function p%. (u>/T). The two spikes correspond 
to the sign of px. changing first from positive to negative, and then positive again. 



(iii) + 2 / dq 



2u 



Aq 
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dr 



n. 



n„ - n„_" ■ 



(l+n q - r )(n q -n r+ u) 



n r n w 



qr 
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Final result 

In analogy with the other master integrals, the only thing left is to collect our final result 
for the function p x . Adding up all the analytic parts and performing the numerics for the 
remaining integrals, we get 



(4.7r) 3V 2 



23 46, A 2 41 . ;mN 
+ — In — + — -px.(u/T) 



18e 18 oj 2 



27 



while the behavior of px. is displayed in fig. 8. 
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